Inelastic Tunneling Spectroscopy in a D-wave Superconductor. 



A.V. Balatsky, Ar. Abanov, and Jian-Xin Zhu 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 

(Dated: February 1, 2008) 

We propose a mechanism to use inelastic tunneling spectrosopy STM to detect a single spin in a 
d-wave superconductor and in a pseudogap state, based on a direct exchange coupling J between 
the surface electrons and the local spin S in a magnetic field. This coupling will produce a kink in a 
dl jdV characteristic at Zeeman energy of the spin Uo- We find that for relevant values of parameters 
signal scales as dl 2 /dV 2 ~ ( JNo) 2 Q(eV — u>o) and could be in the range of 10 -2 of the bare density 
of states where No is the density of states for surface electrons. Scattering in superconductor with 
the coherence peak at gap maximum A leads also to strong features at A + ujo- This suggests a 
new technique for a detection of a local spin excitation with STM. We also consider a detection of 
a local vibrational mode as a simple extension of the spin case. 
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Inelastic electron tunneling STM spectroscopy (IETS- 
STM) is a well established technique that has been 
proven, starting with important experiments of Stipe et 
al. 0. In these experiments a step-like feature in tun- 
neling current and local density of states have been ob- 
served. The physical explanation of the effect is straight- 
forward: once energy of tunneling electrons exceeds the 
energy required to excite local vibrational mode, there 
is a new scattering process that contributes to the scat- 
tering of electrons due to inelastic excitation of the local 
mode Q. Similarly, in case of a single impurity spin S 
in external field the localized spin state will be split with 
the Zeeman gap Hluq = 9^bB. If there is a local exchange 
coupling that allows conduction electrons to scatter in- 
clastically off the local spin in an external magnetic field 
then one has a very similar situation of inelastic scatter- 
ing off localized spin levels. Up to date the single spin 
inelastic tunneling spectroscopy in metals with STM has 
not been observed however. One reason that is often 
mentioned is a Kondo screening of a magnetic spin by 
conduction electrons that leads to a more complicated 
response 0- 

Potential applications of the techniques that are sen- 
sitive to the single spin dynamics include studies of a 
single spin Kondo problem 0, studies of magnetic im- 
purity states in superconductors [f| and single spin de- 
tection and manipulation in a context of a solid state 
quantum computing schemes 0, Q ■ Use of STM for spin 
detection is a promising approach that would allow one 
to combine a sub- Angstrom spatial resolution with the 
high electronic sensitivity. 

We propose to study the inelastic electron tunneling 
spectroscopy of a single spin and of a localized vibrational 
mode with STM in a d-wave superconductor. In case 
of spin there is a crucial difference between metal and 
d-wave superconductor that might make the single spin 
observation in IETS with STM more feasible: vanishing 
DOS in d-wave superconductor drastically suppresses or 
even makes the Kondo temperature Tk = for a single 



impurity spin. This leaves one with a simpler problem of 
a scattering off the single unscreened spin in a media with 
linearly vanishing density of states (DOS) N(lu) ~ u>. 

The proposed approach to identify magnetic sites is 
based on the fact that the local density of states in the 
vicinity of a single impurity spin will have a kink-like 
feature with threshold at eV = tkjQ, luq is a Larmor fre- 
quency of a spin. This would be a natural, albeit not 
tried yet in correlated electron systems, extension of the 
single molecule vibrational spectroscopy and could allow 
a single spin detection. Our results can be summarized as 
follows: i) we find that spin produces the kink- like singu- 
larity in the density of states that as a function of position 
with respect to moment site at low T <C luq,lu <C A is : 



SN(r,w)/N ^2n 2 (% ln(^)) 2 
-lu )(N JS) 2 A(v) , 
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where A(r) is a combination of Greens function of an 
electron in real space (assumed 2D), describing Fricdcl 
oscillation clS cl function of r and reflecting four-fold 
anisotropy due to superconducting gap. The singularity 
at lu = luq(B) changes as a function of applied external 
magnetic field since lu = g^LgB, fig being the Bohr mag- 
neton and g the gyromagnctic ratio; We also find a strong 
feature at energies A + luq as a result of sharp coherence 
peak in DOS of a superconductor, lu ~ —A: 
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(see Eq. 0); these kink-like singularities in dl/dV char- 
acteristic of a STM tunneling in a vicinity of magnetic 
site lead to a step in d 2 I/dV 2 ; ii) the strength of the ef- 
fect is of second order in a dimensionless coupling N J. 
If we take typical values of J ~ 1 — 0. leV and N = 1/eV 
we find that the magnitude of the correction to DOS is 
on the order of 10 -2 . iii)The proposed effect can also 
be trivially expanded to be tried for a molecular spec- 
troscopy of a local vibrational mode. One would have 
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to assume uiq be the eigenfrequency of a local mode, in- 
dependent of the field, and replace JS by the coupling 
constant to the local vibrational mode in Eq. QJ. 

Important potential application of the proposed tech- 
nique is the study of the induced magnetic moment near 
Zn and Ni impurities in the high-T c materials. It has 
been argued that the Zn, Ni and Li impurities in the 
Cu-0 plane generate uncompensated spin. A Zn 2+ that 
substitutes Cu 2+ has a closed shell and is nonmagnetic 
and induced moment near Zn has to be a collective re- 
sponse of the neighbor sites. Claims also have been made 
about evidence of the Kondo effect ||. Direct and in- 
dependent test of magnetism induced by Zn, Ni, and 
Li impurities would be important for our understanding 
of the physics of strong correlations in Cu-0 planes in 
high-T c compounds. The IETS STM technique would 
allow the direct and alternative approach to distinguish 
between magnetic and nonmagnetic sites in a d-wave su- 
perconductor. 

Assume that we have localized magnetic atom with 
spin S on a surface of a d-wave superconductor. Electrons 
in a superconductor interact with the localized spin via 
point-like exchange interaction at one site JS • or: 

H = J2 c U(k)c kCT + ^(A(k)4 T ct_ ki + h.c.) 

k k 

c k<r' T ^' c k'cr' + .9^bS • B , (3) 

k,k',a\<r' 

where Ckc is annihilation operator for the conduction 
electron of spin a, e(k) is the energy of the electrons, 
A(k) = A (cos k x — cos k y ) is the d-wave superconducting 
gap of magnitude A ~ 30meV in typical high-T c mate- 
rials. The local spin S is a \S\ = 1/2. Wc focus here 
on the effect of the Zeeman splitting of the otherwise de- 
generate local spin state in the external magnetic field B 
with splitting energy ujq = lul = g/isB. Below we use 
a mean field description of superconducting state at low 
temperatures T <C T c . Assuming field B <C H C 2 we will 
ignore the orbital and Zeeman effect of the field on the 
conduction electrons H|. 

We are interested in a local effect of inelastic scattering 
of electrons. Thus only local properties will determine 
the conduction electron self-energy. Results we obtain 
will also hold for a normal state with linearly vanishing 
DOS, such as a pseudogap state of high-Tc superconduc- 
tors. In the case of a normal state one would model nor- 
mal pseudogap state with a single particle Hamiltoninan 
H o = Ek c L e ( k ) c k<r with N(w) ~ uj. 

Because of the vanishing DOS in a d-wave supercon- 
ducting state Kondo singlet formation occurs o nly for a 
coupling constant exceeding some critical value [lOj . For 
a particle-hole symmetric spectrum Kondo singlet is not 
formed for arbitrarily large values of J. Another situa- 
tion where Kondo effect is irrelevant is the case of ferro- 
magnetic coupling J. This allows us, quite generally, to 



consider a single spin in a d-wave superconductor that is 
not screened and we ignore the Kondo effect. 

In the presence of magnetic field B||z spin degener- 
acy is lifted and components of the spin S||z and S _L B 
will have different propagators. It is obvious that only 
transverse components of the spin will contain informa- 
tion about level splitting at luq = lul ■ We have therefore 
focused on S + , S~ components only. The propagator in 
imaginary time r is x( T ) = {TtS + (t)S~ (0)} with Fourier 
transform and continuing to real frequency Xo(k>) = 

^-(u+iS) 2 • ^ 0r ^ ree s P' n we nave = tanh(cjo/2T)/2. 
For more general case of magnetic anisotropy this does 
not have to be the case. To be general wc will keep (S z ). 

We begin with evaluation of the DOS correction due 
to coupling to localized spin. Self-energy correction is: 

E(w,) = J 2 t G ( k > - n n ) x + ~ (a,) , (4) 

where G°(k,w,) = [iw, - e(k)][(^) 2 - e 2 (k) - A 2 (k)]- 1 
is the particle Green's function in d-wave superconduc- 
tor, G- 1 = GW- 1 - E, F°(k, w; ) = [A(k)][(^,) 2 - 
e 2 (k) - A 2 (k)]" 1 ; = 2nlT is the bosonic Matsub- 
ara frequency and u>i = (21 + l)wT;l = 0,1,2... is the 
fcrmionic frequency. Using spectral representation and 
analytical continuation onto real axis icu n uj + iS wc 
find for imaginary part of self energy T,(uj) : 

ImE(ijj) = —J 2 (S z )ImG(Ld~LUo){nF(Lu—Ldo)—nB{uo) — ]-} , 

(5) 

where n F (uj) = l/[l + cxp(/3w)], n B (ui) = l/[exp(/3w) — 1] 
are Fermi and Bose distribution functions. This lo- 
cal self-energy leads to the modifications of the DOS. 
In this solution we treat the self-energy effects in G 
to all orders, i.e. G in Eq.lJSJ is full Green's function 
G _1 = Gq 1 — and solution for S is found self- 

consistently for a local vibrational mode. The mod- 
ifications of the superconducting order parameter and 
bosonic propagator were ignored in this calculation. Re- 
sults are presented in Fig. ^ To proceed with analytic 
treatment, unless stated otherwise, we limit ourselves be- 
low to second order scattering in S. Difference between 
self-consistent solution and second order calculation are 
only quantitative and small for small coupling. Correc- 
tions to the Green's function G(r,r',w) = G°(r, r',w) + 
G°(r, 0, w)E(w)G°(0, r', w) + F°(r, 0, w)E(w) F*°(0, r, lu). 
For simplicity we define K(T,uj,ljq) = — [np(w — u>o) ~ 
nsi^o) — 1 — ©(ijJ — wo)] which becomes a step function 
at low T too, the limit wc will focus on hereafter. Cor- 
rection to the local density of states as a function of posi- 
tion comes from the correction to the bare Green's func- 
tion G° : 6N(t,w) = l/7rIm[G°(r,0,cj)I](w)G (0,r,cj) ± 
F°(r, 0, w)E(w)F*°(0, r, lu)], where keeping it general, 
the plus sign corresponds to the coupling to the local 
vibrational mode and minus - to the spin scattering re- 
spectively. The strongest effect will be at the impurity 
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site. For on-site density of states we have: 



SN{r = 0,uj) _ 7T 



2 W — tdQ 



N 



Y (JSN y A 



K(T,(j,(j Q ) 



(6) 



SN(r = 0,w) = 2?r 2 (J5iVo)2ir(r)W)Wo)ln2 / I- - Al 



N 



4A 



x In 



4A 



led + W - A| 



+ (w -> -wo), W ~ |A| , 



(7) 







1 ' ra =°-2 : 

N(o>) 

dN(o>)/do> — 

^ — A-«» ~ 












co =0_4 _ 
A+QJ 


"«n-TT| i i i, i 




Vl i 


°°o J 


III 1 


co =0.6 _ 


1 1 , 1 1 , 


1 





where we used for on-site Green's function G°(0,0,ui) 
AMttM^t) + """ir)' for td < A and we retained 



only dominant real part of G°. In opposite limit 
cd ~ A we retained only dominant imaginary part of 
G°(0,0,cd) = wrAT(w) = -2iAT ln(J^i). At r = we 
have F°(0,0,w) = 0. Complete DOS N(to) and deriva- 



tive 



dN(m) 
dej 



are shown on Fig. ^ For arbitrary position 
N(r, td) we would have to add a Friedel oscillation factor 
A(r) = [|G°(r,w)| 2 ± |F°(r, W )| 2 ] ~ that 
describes the real space dependence of the Green's func- 
tion on distance for small w < A. Here rjJIkpj. is the 
component of r = (rj_,rii) that is along the Fermi sur- 
face near the nodal point of the gap and r||||kp|| is the 
component perpendicular to the Fermi surface at nodal 
point. Existence of the nodes in d-wave case results in 
the power law decay of A(r) in all directions and it has a 
four fold modulation due to gap anisotropy (See detailed 
discussions in PRB '97 reference in ). The final result 
is our Eqs. 1)1121) . 

It follows immediately that 



S w/w ~ SN(v = 0, V)/N ~ (JSN ) 



2 V- 



^e(y-cdo) 



u dV? 



(JSN ) 2 Q(V -wo) 



FIG. 1: The solution for the DOS (black line) and its en- 
ergy derivative (red line) are presented for a local boson 
mode scattering in a d-wave superconductor. The normal 
self-energy was treated self-consistently as a full solution of 
the Eq. JjJ , ignoring vertex corrections and gap modifications. 
Apart from the feature at lo = loo we get also strong satellite 
peaks at A + ujq that are a consequence of a coherence peak 
in DOS of a d-wave superconductor. These satellites are a 
specific property of a superconducting state and will not be 
present in a pseudogap state. These features are best seen in 
Energy scale is given in units of A, the dimensioneless 
coupling constant is taken to be 1. For comparison we plot 
the results for local mode frequency wo/A = 0.2,0.4,0.6 in 
the first three panels. The lower panel gives the results for 
the asymptotic analytic solution, that assumes too -C A using 
Eq.(l), for wo = 0.4. The overall features are similar for both 
cases, however the analytic solution shows a somewhat larger 
feature. 



N{lo) = l/7rImG°(0,0,cd) = (cd/A) 7 Ao with power 7 > 
that is determined by the microscopic properties of the 
(8)material. Then, from Eqs. (I'iKill we have for w C A: 



Here we have used the fact that the derivative of (td — 
cdo)6(cd — u>o) with respect to td yields Q(u — u>o). Thus in 
a d-wave superconductor and in a metal with vanishing 
DOS AT(cd) = Nq^ one should expect a step discontinu- 
ity in d 2 I/dV 2 at the energy of a local mode with the 
strength J 2 Nq (see Fig. This result is qualitatively 
different from the case of conventional metal. For metal 
with energy independent DOS we have from Eq. 10 for 
T<td 
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dV ' dV 



SN(r 



0,V)/N o 
w ) 



(v - w )Te(v 



^-(v-cdor-w 



w ) 



w ) , 
(10) 



dl 
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SN(r = 0, V) ~ J 2 A 3 9(V - w ) 



(9) 



and the second derivative will reveal a delta function 
d 2 I/dV 2 ~ J 2 N$5(lu - lvq) The effect in d-wave super- 
conductor is clearly smaller than correction to DOS in a 
normal metal with the same coupling strength. 

For completeness we also have calculated the effect of 
inelastic scattering in a metal with the more general DOS 



Depending on the value, we get divergent singularity at 
Wo for 7 < I, or a power law rise for 7 > 1. In case of 
7 = 1 we recover the result for d-wave superconductor 
and for a pseudogap normal state. 

Quite generally one can express the results in terms of 
the spectrum of superconductor. We can write ImS(w) 
using spectral representation for G(r, w). In supercon- 
ducting case, using Bogoliubov u a (r) 7 v a (r) for eigenstate 

a, we have G(r,w) = ZJ^0r s + T ^ing 
imaginary part of G(r, td) we arrive for T <C tdo at: 

ImS(w) = ^<S z )[Mr = 0)| 2 <S(td-td -£ Q ) 

^Wo 

+K(r = 0)| 2 <5(w-w + £ Q )] ,w>0. (11) 
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FIG. 2: Appearance of the satellite peaks for an impurity 
resonance u>i mp at u>i mp ± ujo is shown schematically. The 
satellites will have different spectral weight. Imagine we inject 
into system an electron at energy Ui mp +ujq. To create a peak 
at coimp one need to excite local mode and the energy of the 
electron will be equal to the difference between local state 
and local mode energies. Similarly, to obtain the peak at 
Wi m p from injected electron at energy u>im P — IjJo one needs 
to add local mode energy to an electron. For this process 
to occur the local mode has to be excited to begin with and 
hence this process will have very low weight at low T. These 
two processes will also have different matrix elements. Overall 
relative weight of the side peaks is proportional to J 2 Nq which 
we assumed to be small. In case of magnetic scattering, when 
wo = the splitting will be tunable by the field. 



At negative w < one has to replace luq — * — wo in 
Eq. (|llfl . For example, consider a magnetic impurity res- 
onance in d-wave superconductor at energy uJi mp , such 



[Tllllll- Then only the term 



as a Ni induced resonance 
with resonance level E a = Ei mp will dominate the sum 
over eigcnstatcs a in the vicinity of impurity site. In- 
clastic scattering off this impurity induced resonance will 
produce additional satellite split away from the impurity 
level by wo, see Fig- EI Sharp coherence peaks will also 
produce split satellites. Again, for a local phonon mode 
one gets a similar splitting of impurity level with wo now 
being the phonon energy. 

Our results suggest the possibility of single spin detec- 
tion as one monitors the feature in d 2 I /dV 2 as a function 
of position and external magnetic field. If we take exper- 
iment allyseen DOS N ~ 1/eV with JN ~ 0.14, A = 
SOmeV hj and assuming the field of ~ 10T we have 
luq = lmeV (corresponding to the Zeeman splitting of 
~ lmeV in a magnetic field, we have from Eqs. 10181) 



5N(r = 0,w)/Y o ~ 10" 
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We point out here that result is expressed in terms of the 
relative change of DOS of a metal No. For observation of 
this effect one would have to sample DOS in the vicinity 
of eV = ujo oc B. Assuming w — wo = ujq we have from 
Eq. (fl"2*|) S-^/^y ~ 10~ 2 . Expressed as a relative change 



of DOS of a superconductor N(uj) 
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N u>/A effect is: 

2 ^^e(u! - w ). 



brational modes of localized molecules in inelastic elec- 
tron tunneling spectroscopy STM, IETS-STM Q. The 
satellites at A + w produce the effect on the scale of 
unity and clearly seen even for small coupling. The im- 
portant difference is that for localized spin the kink in 
DOS is tunable with magnetic field and this should make 
its detection easier. 

In conclusion, we propose the extension of the inelastic 
tunneling spectroscopy on the strongly correlated elec- 
trons states, such as a d-wave superconductor and pseu- 
dogap normal state. The DOS in these systems has a 
nontrivial energy dependence of general form N(lo) ~ 
w 7 , 7 > 0. This technique could allow for a Zeeman level 
spectroscopy of a single magnetic center, thus, in princi- 
ple, allowing a single spin detection. We find the feature 
in dl/dV ~ (w — wo) 7_1 0(w — uiq) near the threshold 
energy wo- We also find strong satellite features near the 
gap edge due to coherence peak for a superconducting 
case. The singularity is a power law and qualitatively dif- 
ferent from the results for a simple metallic DOS 0. For 
the relevant values of parameters for high-Tc the feature 
is on the order of several percents and makes the fea- 
ture observable in these materials. Similar predictions 
are also applicable to the local vibrational modes, where 
wo becomes a vibrational mode frequency. 
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